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ABSTRACT 


Results  of  the  calculation  of  helicopter  airloads  using  a  lifting  sur¬ 
face  theory  solution  are  presented  and  compared  with  experimental  data.  These 
results  indicate  that  a  very  accurate  wake  geometry  model  will  be  required  in 
order  to  make  full  use  of  the  accuracy  of  the  lifting  surface  theory  solution. 
It  is  observed  that  the  experimental  vortex  induced  loads  decrease  as  the 
vortex  moves  inboard  along  the  blade;  this  phenomenon  in  the  behavior  of  the 
vortex  wake  of  a  rotor  requires  more  Investigation. 
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Rotor  thrust  coefficient: 


tiirust 
P(nR)  ^TTR^ 


Rotor  radial  coordinate 

Blade  radius 

Root  collective  pitch 

Advance  ratio:  forward  speed  divided  by  rotor  tip  speed 
Air  density 


Solidity  ratio: 


total  blade  area 


TTR 


V  Rotor  aziauth  coordinate 

Q  Rotor  rotational  speed 


INTRODUCTION 


A  lifting  surface  theory  solution  has  been  developed  for  a  model  prob¬ 
lem  for  vortex  induced  airloads  (reference  1) .  The  model  consists  of  an  in¬ 
finite  aspect  ratio  wing  in  a  subsonic,  compressible  free  streaun,  and  a 
straight,  infinite  vortex  at  an  arbitrary  angle  with  the  wing  (figure  1) .  Using 
the  exact  numerical  results  from  linear  lifting  surface  theory,  an  approximate, 
closed  form  solution  for  the  vortex  induced  loads  was  obtained.  The  solution 
in  this  form  is  suitable  for  application  to  the  calculation  of  rotary  wing  air¬ 
loads)  procedures  for  this  application  were  developed  in  reference  1.  This 
report  presents  results  of  calculations  using  the  lifting  surface  theory  solu¬ 
tion  for  the  vortex  induced  loads,  and  compares  the  theoretical  calculations 
with  available  experimental  data. 


THE  CALCULATION  OF  HELICOPTER  AIRLOADS 

Exoerimental  data  is  available  (references  2  and  3)  from  flight  tests 
of  a  four-blade  rotor  of  a  Sikorsky  H-34  helicopter;  a  description  of  this 
rotor  and  its  instrumentation  may  be  found  in  reference  2.  Five  cases  were 
chosen  for  study  (flights  number  7,  12,  17,  and  25  from  reference  2,  and  the 
U  ■  0.18  case  of  reference  3);  these  cases  were  chosen  because  they  show  vor¬ 
tex  induced  loads  at  r  ■  .95. 

Figures  2  through  6  present  a  comparison  of  calculated  and  experimental 
section  lift,  for  radial  stations  r  ■  .95,  .85,  .75,  .55,  and  .25,  respectively. 
The  experimental  data  is  from  flight  number  7  of  reference  2;  this  case  has  an 
advance  ratio  of  M  ■  0.15.  The  calculations  used  the  lifting  surface  theory 
solution;  a  rigid  wake  geometry  was  used  and  the  values  of  the  wake  inflow  and 
the  first  harmonic  flapping  were  adjusted  to  give  agreement  of  the  loads  at 
the  r  ■  .95  station.  The  inboard  loads  tend  to  be  low;  the  calculated  tiurust 
gives  C^O  ■  0.082,  while  the  experimental  value  is  about  C^o  «  0.089.  The 
greatest  discrepancies  in  the  prediction  of  the  inboard  loads  arise  because 
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of  the  use  of  a  rigid  wako  geometry.  Figures  2  and  6  also  include  the  re¬ 
sults  of  a  calculation  using  a  nonrigid  wake  geometry,  with  a  large  vortex 
core  radius  to  account  for  local  distortion  due  to  the  vortex/blade  inter¬ 
action.  The  non-rigid  wake  geometry  was  supplied  by  M.P.  Scully,  using  a 
development  of  his  method  described  in  reference  4  for  the  calculation  of  the 
self-induced  distortion  of  the  tip  vortices  In  the  rotor  wake.  The  prediction 
of  the  inboard  loads  is  substantially  improved  by  the  use  of  the  nonriqld  wake 
geometry.  Figures  7  through  9  present  polar  diagrams  of  the  section  lift,  for 
the  experimental  data,  the  rigid  wake  calculation,  and  the  nonrigid  wake  cal¬ 
culation,  respectively.  It  is  seen  that  a  good  wake  geometry  model  Is  neces¬ 
sary  in  order  to  predict  the  rotor  loads.  While  the  lifting  surface  solution 
allows  an  accurate  calculation  of  the  loads.  Its  use  requires  the  distribution 
of  the  downwash  over  the  rotor  disk,  and  the  downwash  must  be  calculated  from 
the  vortex  wake  of  the  rotor.  Moreover,  a  verv  high  degree  of  accuracy  Is 
required  in  order  to  correctly  predict  the  closeness  of  the  vortex  to  the 
blade,  to  which  the  downwash  Is  most  sensitive. 

Figures  10  through  12  present  a  comparison  of  calculated  and  experi¬ 
mental  section  lift,  for  radial  stations  r  ■  .95,  .85,  and  .75,  respectively. 

The  experimental  data  is  from  the  U  ■  0.18  case  of  reference  3.  The  calcula¬ 
tions  used  the  lifting  surface  theory  solution;  a  rigid  tmke  geometry  was  used 
and  the  value  of  the  wake  Inflow  adjusted  to  give  a  correct  level  of  the  peak- 
to-peak  vortex  Induced  loads  at  r  *  .95.  Flapping  was  calculated  theoretically 
(rather  than  adjusted  to  give  agreement  of  the  loads  at  r  ■  .95),  producing 
some  discrepancy  In  the  loads.  In  order  that  the  theoretical  peak-to-peak 
vortex  loads  on  the  advancing  side  of  the  disk  (which  are  due  to  the  tip  vor¬ 
tex  from  the  preceding  blade)  decrease  as  the  experimental  data  show.  It  was 
necessary  to  arbitrarily  push  t)ie  vortex  fartlier  from  the  blade  as  the  blade 
passed  over  It. 

The  use  of  the  lifting  surface  theory  solution  does  require  more  calcu¬ 
lation  than  the  use  of  lifting  line  theory,  but  still  the  airloads  calculation 
remains  small  compared  with  the  downwash  calculation.  The  times  to  calculate 
the  down%#ash  at  one  point  on  the  rotor  disk  due  to  one  line  element  In  the  wake, 
to  locate  the  point  of  nearest  approach  of  a  tip  vortex  to  the  blade,  and  to 
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calculate  the  load  at  one  point  on  the  disk  due  to  one  vortex  using  the  lift¬ 
ing  surface  solution  are  in  the  ratio  1:28:22.  In  a  typical  case,  for  one 
cycle  of  the  blade  around  the  disk  the  total  tines  to  calculate  the  downwash, 
to  locate  the  points  of  nearest  approach  of  the  vortices,  and  to  calculate  the 
loads  using  the  lifting  surface  solution  were  in  the  ratio  25:1:3. 

The  feature  of  the  airloads  referred  to  above  as  local  distortion  due 
to  vortex/blade  interaction,  requiring  either  a  large  core  size  or  pushing  the 
vortex  away  to  obtain  correct  inboard  loads  due  to  the  vortex,  is  actually 
much  more  involved.  Figures  13  through  17  compare  the  experimental  section 
loads  for  the  five  cases  studied.  It  is  seen  that  for  every  case  the  effect 
of  the  vortex  is  greatly  reduced  as  it  moves  inboard  along  the  blade.  Current 
wake  geometry  models  indicate  that  on  the  advancing  side  of  the  disk  the  tip 
vortex  of  a  blade  remains  very  close  to  the  plane  of  the  rotor,  and  is  pushed 
dotmward  after  the  passage  of  the  following  blade.  The  vortex/blade  separa¬ 
tion  is  small  and  varies  little  as  the  blade  moves  over  the  vortex.  At  the 
advance  ratios  of  the  cases  considered  here  the  vortex  moves  inboard  at  least 
to  the  r  >  .70  station.  With  this  geometry  the  vortex  induced  loads  should 
rmsain  about  the  same  as  the  vortex  moves  inboard;  instead  the  loads  are  ob¬ 
served  to  decrease  substantially.  The  present  calculation  procedure  might 
account  for  this  effect  by  having  the  blade  push  the  vortex  away  as  it  passes 
over  it,  or  by  using  a  large  viscous  core  size.  While  these  two  features  are 
undoubtedly  a  part  of  the  phenomenon,  the  calculations  above  indicate  that  the 
vortex  would  have  to  be  pushed  too  far  away,  or  too  large  a  core  size  would 
have  to  be  useu  for  either  of  these  to  be  the  entire  cause.  Other  possible 
causes  are  vortex  bursting  due  to  the  presence  of  the  blade,  or  the  inter¬ 
ference  of  the  vortex  with  the  wake  vorticity  it  generates  behind  the  blade 
propagating  up  the  vortex.  For  a  very  close  vortex  the  viscous  core  will  be 
in  contact  with  the  boundary  layer  of  the  blade;  moreover,  large  radial  ve¬ 
locities  will  be  induced  on  the  blade  by  the  vortex;  these  are  likely  important 
features  of  the  vortex  behavior.  The  nature  of  this  phenomenon  and  its  causes 
are  at  present  unknoim.  Zt  is  evident,  hoiiwver,  that  it  not  sufficient  to 
consider  the  rotor  wake  simply  as  composed  of  well-behaved  vortex  lines  and 
sheets.  There  is  yet  a  great  deal  to  be  learned  about  the  nature  of  the  wake 
of  a  rotary  wing. 
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FIGURE  I  MODEL  PROBLEM  FOR  VORTEX 
INDUCED  AIRLOADS 
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FIGURE  2  COMPARISON  OF  SECTION  LIFT  CALCULATED 
USING  LIFTING  SURFACE  THEORY  WITH 
EXPERIMENTAL  RESULTS  (FLIGHT  No.  7 
OF  REF  2)‘  ^«QI5,  ^®«I2.5®  ,  r  *,95 
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FIGURE  3  r  «  .05 
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FIGURE  4  r  »  .75 
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FIGURE  5  r  »  .55 


9 


1 


10 


180 


FIGURE  8  LIFTING  SURFACE  THEORY  RESULTS  WITH 
RIGID  WAKE  (LB/IN) 
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SECTION  LIFT  (  LB/IN) 


FIGURE  10  COMPARISON  OF  SECTION  LIFT  CALCULATED 
USING  LIFTING  SURFACE  THEORY  WITH 
EXPERIMENTAL  RESULTS  (FROM  REF.  3)« 
ft  *0.18,  8*»l3i7*  .  r  *.95 
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SECTION  LIFT  (LB/IN) 
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SECTION  LIFT  (LB/IN) 


FIGURE  I2«  r*.75 
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FIGURE  13  EXPERIMENTAL  SECTION  LIFT 
(FLIGHT  No.  7  OF  REF.  2)  ‘ 
/i*O.I5.  ^•■12.5®,  Cj/o-s  0.089 
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FIGURE  14  EXPERIMENTAL  SECTION  LIFT 
(FLIGHT  No.  12  OF  REF  2  )« 
fi,  «0.25.  8**  13.4*,  Cj/o-  «  0.091 
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FIGURE  15  EXPERIMENTAL  SECTION  LIFT 
(FLIGHT  No.  17  OF  REF.  2)  • 
fi  *  0.26,  12.7®,  Cj/a  =0.068 
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FIGURE  16  EXPERIMENTAL  SECTION  LIFT 
(FLIGHT  No.  25  OF  REF.  2)  « 

/I  «  0.08,  12.7*.  Cj/a  ■  0.072 
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SECTION  LIFT  (LB/IN) 
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FIGURE  17  EXPERIMENTAL  SECTION  LIFT 
(FROM  REF  3)* 

*  0.18,  8®*  13.7®,  Cy/er  =  0.089 
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